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ABSTRACT:  The  upward  ndgration  of  pulsating  explosion  bubbles  under  the 
Influence  of  buoyancy  (gravity)  cannot  be  studied  on  a  model  scale  in  free 
water,  but  only  in  test  tanks.  If  a  stationary  test  tank  is  employed  the 
air  pressure  above  the  water  surface  must  be  reduced  to  assure  similitude. 
This  report  deals  with  another  alternative,  namely  a  tank  which  is  subjected 
to  a  suitably  selected  acceleration  during  the  eiiploslon  test.  The  case 
considered  is  that  where  the  acceleration  is  produced  by  a  centrifuge. 

The  Coriolis  accelerations  distort  the  motions  of  the  migrating  bubble  in 
two  ways.  First,  the  "upward"  (i.e.  towaz^  the  water  sxirface)  motion  of  the 
bubble  is  deflected  to  the  side.  This  sideward  motion  can  be  observed 
experimentally.  In  the  second  place,  the  rotation  affects  the  magnlttide 
I  of  the  "upward"  motion  and  introduces  a  systematic  error  which  cannot  be 

detected  by  e:^rlmental  methods. 

The  bubble  migration  in  a  rotating  system  is  treated  using  the  simplifying 
assus^tlon  that  the  bubble  retains  its  spherical  shape.  The  "upward"  as 
well  as  sideward  migration  is  calculated  for  the  first  and  second  cycle  of 
the  pulsation.  In  a  typical  example,  the  systematic  error  in  the  "upward" 
migration  was  about  1^,  which  is  probably  less  than  the  e:]q>erimental  accuracy 
of  such  studies. 
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NmRATION  OF  EXPLOSION  BUSBIES  IN  A  BOTASTUK}  TEST  TANK 

The  design  of  a  high  gravity  tank  for  the  study  of  underwater  e^loslon 
phenomena  Is  greatly  slsqpllfled  If  the  tank  Is  accelerated  by  a  centrifuge 
Instead  of  a  linear  accelerator.  This  report  shows  that  the  rotation  does 
not  Introduce  errors  which  exceed  the  e^^rlnental  accuracy  of  such  tests. 
Therefore,  the  rotating  tank  Is  a  promising  tool  for  model  experiments  on 
migrating  explosion  bubbles. 

This  paper  is  the  outcome  of  calculations  which  were  made  during  the  planning 
studies  for  a  rotating  high  gravity  tank.  Since  the  background  as  well  as 
the  details  of  these  calculations  are  of  Interest  to  the  users  of  such  tanks, 
they  are  made  available  In  this  report. 

This  work  was  done  under  Task  NOL-440,  and  was  supported  by  the  Defense 
Atomic  SiqE>Port  Agency. 
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I.  INTRODUCTION 

Various  proposals  have  been  made  in  the  past  to  study  the  upward  motion 
of  underwater  explosion  bubbles  on  a  model  scale  in  a  tank  which  is  subjected 
to  a  hi^  acceleration,  to  By  this  way  a  gravitational  field  can  be 

established  which  assxrres  similitude  of  the  bubble  motions  induced  by 
buoyancy.  To  some  extent  this  objective  can  be  also  fulfilled  in  a  "vacuum 
tanlc".  Reduction  of  the  air  pressure  above  the  water  can  provide  the 
conditions  for  similitude  of  bubble  migration.  This  method  requires  that 
certain  additional  conditions  are  satisfied,  namely  either  a  specific  water 
temperatiire  so  that  a  controlled  boiling  occurs  at  the  sirrface  of  the  bubble 
or  the  use  of  a  liquid  of  low  vapor  pressure  and  special  explosives  of  low 
density.  In  both  cases  it  is  necessary  to  rely  heavily  on  theory  in  order 
to  specify  the  test  conditions  ^ .  Such  a  study  can  hardly  be  considered 
as  a  true  experiment,  but  rather  as  a  combination  of  a  theoretical  and 
e:}p)erimental  approach.  Nevertheless,  quite  satisfactory  results  have  been 
obtained  with  this  method  ^  and  GJ . 

Since  shockwave  and  bubble  pulse  phenomena  are  not  scaled  by  the 
vacuum  tank  technique,  this  method  fails  for  phenomena  where  pressure  effects 
as  well  as  bubble  migration  are  important.  Examples  are:  effect  of  the 
migrating  bubble  on  targets,  interaction  of  the  bubble  with  the  water  surface, 
and  surface  phenomena.  (Studies  of  the  latter  effects  will  require  the  use 
of  special  media  in  the  tank,  for  instance  the  use  of  a  layer  of  dust  at  the 
water  surface  separated  from  the  water  by  a  membrane.  This  would  also  make 
possible  the  study  of  base  surge  phenomena. )  In  a  vacuum  tank  it  is  also 
impossible  to  simulate  the  migration  and  condensation  phenomena  of  steam 
bubbles  such  as  are  produced  by  nuclear  underwater  explosions.  For  all  these 
pinposes  a  hi^  gravity  tank  is  needed,  as  discussed  in  j/. 

The  high  gravity  tanks  originally  proposed  employed  linear  accelerators. 
Subsequently,  R.  S.  Price  (NOL)  has  proposed  consideration  of  a  tank  mounted 
on  a  centrifuge  and  pointed  out  the  great  advantages  in  structural  design, 
ease,  and  safety  of  operation  of  this  arrangement  and  the  practically 
unlimited  duration  available  for  the  explosion  test.  Tests  yielding 
encouraging  results  have  been  miade  with  a  small  centrifuge  at  the  Naval 
Ordnance  Laboratory  _2/  a  large  centrifuge  at  Sandia  Corporation, 

New  Mexico.  Evalxiation  of  the  latter  tests  showed  satisfactory  agreement 
with  the  bubble  migration  theory  presented  in  j/.  Althoi;i^  this  theory 
reproduces  all  existing  experimental  evidence  with  good  accxiracy,  this 
comparison  is  hardly  sufficient  to  demonstrate  that  bubble  migration  can  be 


mU2R  6l-lk^ 


Ir 

properly  siimilated  in  a  rotating  tank.  In  particular,  it  is  not  obvious  that 
motions  resulting  from  the  Coriolis  accelerations  could  not  adversely  affect 
the  migration  in  such  a  tank  and,  thus.  Introduce  a  systematic  error  ■vAiich,  ^ 

on  the  basis  of  an  empirical  comparison,  would  remain  undetected.  It  will 
be  shown  below  that  this  error  is  small  and  that  Price’s  argximent  is  correct 
that  all  detrimental  effects  caused  by  the  rotation  decrease  with  increasing 
length  of  the  arm  on  •vftxich  the  test  tank  is  mounted. 

II.  EQUATIONS  OF  MOTION  IN  A  ROTATING  SYSTEM 

We  assme  a  system  rotating  about  a  vertical  axis  with  constant  speed 
and  consider  the  motions  in  the  plane  of  the  rotation  only.  There  are  two 
coordinate  systems  of  interest  to  us:  the  stationary  coordinate  system  X,Y 
and  the  rotating  one  x,y.  (The  rotating  coordinate  system  is  that  which 
appears  in  the  photographic  pictirres  of  the  rotating  tank  taken  with  a  camera 
fixed  to  the  rotating  system. )  The  X-  and  x-  coordinates  measure  the  distance 
from  the  center  of  the  rotation,  the  positive  direction  being  away  from  the 
center.  The  Y-  and  y-  coordinates  are  perpendicular  to  X  and  x,  positively 
directed  in  the  direction  of  rotation. 

.  ^ 

Although  the  following  equations  can  be  given  in  the  most  concise  form 
by  means  of  vectors,  we  will  write  the  equations  also  in  Cartesian  coordinates. 

The  coordinates  in  the  moving  system  are  given  in  terms  of  those  of  the  ♦ 

resting  system  by 


or 


the  velocity 


or 


X  =  X  cosmt  +  Y  sinut 
y  =  Y  coscut  -  X  sincut 


X*  =  X- 

coscut  + 

Y*  sinciot 

+  ym 

y.  =  Y- 

coscut  - 

X*  sirKut 

-  xa> 

>  ^ 

to  =  JR 

dt  dt 

1  1 

(2) 


A 
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the  acceleration 


X*  •  = 

X*  *  COSCDt 

+  Y** 

sinwt 

+  xiD^  +  2  y'o) 

y  = 

Y*  *  cosoit 

-  X‘* 

sinxit 

2 

+  yu)  -  2x*a) 

= 

dt^ 

.  [■ 

dt^ 

2c^ 

^  • 

dt  - 

The  angular  velocity  is  denoted  by  cu.  (The  vector'^  is  perpendicular  to  the 
x-y  and  X-Y  plane,  but  the  vectors  and  |^^['^'^]  llie  vithin  this  plane. 

Figure  1.  According  to  the  rules  of  the  vector  analysis,  the  x-components  of 
[lu'r]  andi  ^  are  -yco  and  respectively?  the  y-ccmponents  xro  and 

respectively.) 


Equations  (l)  and  (2)  reflect  the  well  known  fact  that  position  and 
velocity  in  the  moving  system  can  be  obtained  by  means  of  a  superposition 
of  the  motion  relative  to  the  resting  system  and  the  motion  due  to  rotation. 
Such  a  simple  superposition  is  not  applicable  to  acceleration.  A  body  moving 
relative  to  the  rotating  system  behaves  as  if  it  were  subjected  not  only  to 
the  absolute  and  the  centrifugal  acceleration  but  to  an  additional  supple¬ 
mentary  acceleration  -2  [q  .Si  1  called  the  Coriolis  acceleration.  This 
s-upplementaxy  acceleration  it^of  particular  interest  in  oxir  problem. 


If  there  are  no  forces  acting  on  a  body  the  absolute  accelerations  X*  * 
and  Y**  vanish.  A  body  free  to  move,  exhibits  motions  within  the  rotating 
system,  as  if  it  were  accelerated  by  x’ •  and  y,  equation  (3),  with 
X“  =  0  and  Y* '  =0.  Relative  to  the  resting  system  this  body  moves 
imaccelerated,  i.e.  with  constant  or  zero  velocity. 


The  equations  of  motion  of  a  body  free 
rotating  coordinates 

X*  *  -  2y*ai  -  xro^  =  0 

y. .  +  2x*a5  -  yo)^  =  0  . 

The  solution  of  these  differential  equations 

X  =  (xq  +  X^  t)  COSOJt  + 

y  =  (yo  +  Yq  t)  cosoit  - 


to  move  are  in  terms  of  the 

(4) ' 

can  be  immediately  found  from 

(5) 

(xo  +  X*  t)  sinmt  . 


(1) 
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The  subscript  o  refers  to  the  position  and  velocity  of  the  body  at  t  =  0.  The 
velocities  relative  to  the  rotating  system  are 

X*  X*  +  a3y_ 

°  °  (6) 


III.  SPHAY  MOTION 

Before  applying  these  equations  to  the  bubble  migration,  ve  consider 
the  simple  example  of  the  motion  of  the  water  spray  ■vftiich  is  thrown  into  the 
air  by  the  undeimter  explosion.  In  the  case  of  a  test  tank  ■vSiich  rotates 
with  great  speed  arovmd  a  vertical  axis,  the  water  surface  is  essentially  a 
section  of  the  surface  of  a  vertical  cylinder.  The  spray  particles  "iftiich 
are  knocked  off  by  the  explosion  leave  the  water  surface  normally  with  the 
velocity  -Xq,  l.e.  they  move  initially  toward  the  center  of  the  rotation, 
y^  is  zero  for  the  particle.  If  we  ignore  air  resistance,  such  a  particle 
moves  unaccelerated  in  the  resting  system  and  equations  (5)  apply.  With  (6), 
we  obtain 


X  =  (x^  -  Xq  t)  cosojt  +  x^mt  sincot 


y  =  X  cot  coscot  -  (x  -  X*  t)  sincot 
o  00 


(7) 


The  radial  distance  of  the  particle  from  the  center  of  rotation  is 
+  y^  =  Xq^o)^^  +  (l-pcot)^  ,  (8) 

■sdiere 

p>0  .  (9) 

^o  “ 

The  mcjment  Wtien  the  particle  falls  back  on  the  water  surface  is  given 
by  the  time  t  for  vhlch  the  radial  distance  beccxnes  equal  to  the  original 
radius  x^.  This  time,  t^,  is 


CO  t„ 

s 


2  p 

Tvp 


The  function  (lO)  has  a  maximum  for  p  =*  1, 
can  attain  is  tmity. 


1.  •  c  • 


(10) 

the  maximum  value  which  co  t„ 


Figure  2  shows  the  trajectories  of  a  particle  as  observed  from  a  point 
\flilch  rotates  with  the  tank.  (The  water  surface  is  drawn  under  the  assumption 
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of  a  very  vide  tank,  for  instance  a  cylindrical  shell  partly  filled  with  water 
•vdiich  rotates  arcnmd  its  axis.)  The  figure  illustrates  the  sideward  motion 
due  to  the  Coriolis  forces.  In  a  linearly  accelerated  tanlc  all  these 
trajectories  wotild  he  straight  "up"  and  would  coincide  with  the  x-axis.  If 
the  acceleration  of  such  a  tank  lasted  long  enou^,  every  particle  would  fall 
hack  to  the  same  point  -vSiere  it  originated,  even  those  which  are  thrown 
upward  ^Tilth  an  exceedingly  hi^  velocity.  In  a  rotating  tank,  such  a 
particle  hits  the  opposite  water  surface. 


The  spray  motion  as  observed  from  a  fixed  point  of  observation  is  less 
coEiplicated.  Price  (IIOL)  has  made  the  following  calculations  during  the 
planning  phase  of  the  centrifugal  test  tank.  From  the  fixed  point  of  observa¬ 
tion,  the  ira.ter  surface  also  appears  to  be  a  cylindrical  sxxrface.  The  trajec¬ 
tory  of  the  particle  is  a  straight  line  A  B,  Figure  3#  with  the  inclination 


tan  = 


^o  ^ 

X* 

o 


(11) 


The  time  ■which  it  takes  for  the  particle  to  travel  from  A  to  B  corresponds 

to  the  previous  t  .  It  is  obtained  from  s  =  vt  -  x  oA  /sin  From 
_  s  s  o  S' 

Figure  3  it  is  imecSj.ately  seen  that 


s 

sin  ©g 


fo _ 

sin  <l> 


and  we  obtain 


CO  t,,  =  sin  0„  , 

or  since  0  =  2  ijc  (Figure  3) 

s  ^ 


(12) 


(13) 


2  tan  ijc 

CO  t  =  2  sin  cost];  =  - 5— 

1  +  tan'^ 


2  p 


(14) 


This  is  the  same  relationship  as  obtained  from  the  equations  -v^iich  apply  to 
the  rotating  system. 


IV.  BUBBLE  MIGRATION 

The  theoretical  study  of  the  strong  migration  of  explosion  bubbles  is 
considerably  conplicated  by  the  change  of  shape  which  the  pulsating  bubble 
undergoes  ^vhen  it  contracts  to  and  re-expands  from  its  mi m'.rmiTn  size.  (Bubbles 
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■which  migrate  slightly  remain  sphericgj.. )  A  rapidly  migrating  hubhle  is 
almost  exactly  spherical  at  its  first  maximum,  but  it  assumes  a  kidney- like 
shape  or  even  that  of  a  torus  at  the  minimum.  This  cheuige  of  shape  effects 
the  migration  of  the  bubble.  Calculations  assuming  a  spherical  bubble  shape 
throu^out  the  pulsation  have  resulted  in  migrations  too  large  by  about  40^. 
Since  it  is  our  purpose  to  obtain  an  estimate  of  the  possible  deviation  of 
the  migration  caused  by  the  rotation  of  the  tank  such  an  over-estimate  will 
not  be  objectionable,  if  the  calculated  deviations  are  small.  This  would 
mean  that  the  actiial  deviations  are  smaller.  No  theory  exists  today  >diich 
describes  migration  and  bubble  shapes  in  a  satisfactory  manner,  and  it  would 
be  an  exceedingly  difficult  undertaking  to  apply  such  a  theory  to  a  rotating 
system.  Assuming  a  spherical  bubble  shape,  the  migration  in  a  rotating  tank 
can  be  approximately  calculated  in  closed  form. 


As  is  to  be  expected,  rotation  introduces  'two  additional  migration  terms 

which  are  due  to  the  Coriolis  acceleration:  One  is  in  the  direction  of  the 
"gravity”  migration,  the  other  is  peipendicular  to  it.  The  first  of  these 
represents  the  systematic  "error"  of  the  migration  in  a  rotating  tank.  An 
additional  error  occurs  because  the  bubble  migrates  into  a  region  of  different 
centrifugal  acceleration.  These  errors  are  a  part  of  the  "gravity"  migration 
to  be  studied.  They  cannot  be  eliminated  by  experimental  methods.  In  con¬ 
trast,  the  side  migration  can  be  readily  detected  and  eliminated  by  observing 
■the  motion  in  a  plane  perpendicular  to  "the  plane  of  rotation.  It  is  the 
purpose  of  this  study  to  estimate  the  magnitude  of  the  migration  errors. 


The  hydrodynamic  equations  of  an  inviscid  and  incompressible  fluid  in  a 
rotating  system  ao'e  (Squire,  reference  (l2)) 

r7“^ 

V  V  *  0 


VQ  VP/p  «  0, 


(15) 


where  v  is  the  vector  of  the  velocity  and  r  that  of  the  coordinates  relative 
to  the  rotating  system.  Qis  the  potential  of  a  body  force  per  xinit  mass,  P, 
the  pressure,  and  p,  the  density  of  the  fluid.  (15)  can  be  readily  deduced 
from  the  common  hydrodynamic  equation  and  from  (4)  considering  that  Sf^dt  + 
(v  is  the  acceleration  of  a  fluid  particle  in  the  rotating  frame.  The 
acceleration  in  the  hydrodynamic  equation  must  be  referred  to  the  resting 
frame.  Relation  (4)  gives  the  terms  ■which  convert  the  acceleration  from  one 
frame  to  the  other. 
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It  is  generally  assumed  that  the  fluid  motion  of  a  migrating  bubble  is 
irrotational  and  that  a  velocity  potential  exists.  Strictly,  this  is  not  true, 
f  but  it  is  an  acceptable  approximation  for  the  major  portion  6f  the  pulsation. 

In  particTilar,  the  motion  of  a  spherical  bubble  in  an  inviscid  fluid  is 
irrotational  in  a  resting  system.  In  a  rotating  tank,  the  fluid  is  at  rest 
with  respect  to  the  moving  frame  before  the  initiation  of  the  bubble  pulsation. 
Thus,  the  motion  starts  from  the  rest  and  for  spherical  bubbles  in  a  rotating 
system  one  may  also  assimie  the  existence  of  a  velocity  potential  cp- 

If  the  flow  is  axi symmetrical  with  respect  to  the  x-axls,  the  Coriolis 
term  in  (15)  can  be  written  as*(2io^^^)/Tl,  lAiere  ^  is  Stoke’ s  stream  function 
and  T],  the  coordinate  perpendicular  to  the  axis  of  sjrmmetry.  This  form  is 
still  not  suitable  for  the  Bemoialli  equation  >Aiich  is  the  only  possibility 
to  arrive  at  a  solution  in  a  sinple  way.  Hence,  we  replace  this  term  approxi¬ 
mately  by-ScoFiJ''*.  The  properties  of  the  function  yjf*  will  be  discussed 
below.  On  this  basis,  (15)  takes  the  form 

2  o 

V(-  /2  -  2co^*  -  u)^(x  +  y^)/2  +  Q  +  P/p)  =  0.  (l5a) 

Integration  yields  the  Bernoulli  equation  and,  thus,  an  expression  for  the 
pressure  P.  The  components  of  an  external  force  'vftiich  may  act  on  the  sphere 
are 

+1 

»  2tt  J  P(A)  cos  9  d(cos  6) 

(16) 

+1 

■  -  2tt  A^  J  P(A)  sin  6  d(sin  0), 

-1 


where  A  is  the  radius  of  the  sphere  and  P(A)  the  pressure  on  its  surface. 
Since  there  are  no  external  forces  on  a  bubble,  F^  and  F  must  be  zero  at  all 
times.  Evaluation  of  these  integrals  will  yield  the  desired  expressions  for 
the  migration  of  the  bubble. 

The  velocity  potential  cp  of  a  si^iere  which  has  a  variable  radius  A(t) 
and  tftiich  moves  with  the  velocities  Ay*  and  Ax*  in  the  y  and  x  direction 
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respectively  can  te  r^resented  ty  that  of  a  variable  source  and  tvo  dipoles. 
Using  polar  coordinates  r,  9,  the  potential  is 

cp(t,  r,  0)  «  •  cos  0  +  ^  ■  sin  0.  (l7a) 

For  the  corresponding  function  >1^*,  we  set  tentatively 

3  3 

^*(t,  r,  0)  "xjr*  +  sin  0  •  ^  cos  0  .  (iTb) 

s  2r^  2r® 

Since  the  source  does  not  contribute  to  the  migiation,  i|if*  is  not  needed  for 
our  purpose.  ® 

following  properties: 

ax 

(Ifc) 

“  ^  • 


The  function^*  must  have  the 


^  " 

ax  “ 


f 


n 


For  a  plane  flow,  these  requirements  sure  satisfied  by  the  stream  function. 

For  a  three  dimensional  flow,  as  it  is  of  interest  here,  a  function 
probably  does  not  exist.  In  fact,  (l7a)  and  (l7b)  do  not  satisfy  (l7c). 

Still,  (l7b)  is  useful  as  an  approximation.  It  is  readily  verified  that 
the  second  terms  of  (l7a)  and  (l7b)  satisfy  the  first  row  of  (l7c),  but  not 
the  second  row.  Also,  the  third  terms  of  (l7a)  and  (l7b)  satisfy  the  second 
row  of  (rfc),  but  not  the  first  one.  Coraparison  with  (l5)  shows  that  the 
components  of  ■»ftiich  satisfy  (lie)  are  just  those  idiich  determine  the 

Coriolis  acceleration.  Of  course,  this  does  not  mean  that  (l7b)  is  anything 
better  than  an  approximation.  When  inserted  into  the  integrals  (l6)  those 
terms  of  (iTb)  idiich  contribute  at -all  axe  correct  at  the  upper  and  at  the 
lower  limit,  but  faulty  at  the  mid-point  of  the  range  of  integration.  However, 
at  this  point  the  integrant  vanishes  because  of  the  sin  0  or  cos  0  under  the 
Integred.  Thus,  the  error  is  zero  at  both  limits  and  at  the  mid-points  of 
the  integie^tiott.  Since  the  integrants  of  (l6)  contribute  most  strongly  near 
the  limits,  Where  is  most  nearly  correct,  the  use  of  (l7b)  is  justified 
as  an  approximation. 
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In  the  X  -  y  coordinate  system,  used  so  far,  the  positive  x-axis  points 
avay  from  the  center  of  the  rotation.  The  water  in  a  tank  fixed  to  this  axis 
will  he  pressed  in  the  positive  x-direction  and  the  "gravity"  migration  will 
he  in  the  negative  x-direction.  In  underwater  explosion  work,  the  z- coordinate 
commonly  refers  to  the  vertical  coordinate  pointing  toward  the  water  surface. 
The  origin  of  z  is  at  the  point  of  explosion.  Denoting  the  "xqpward"  migration 
hy  Az,  we  have 


X  ■  R  -  Az  +  r  cos  0, 

where  R  is  the  di&tsmce  between  the  center  of  rotation  and  the  i)oint  of 
explosion.  The  y-coordinate  points  in  the  "horizontal"  direction,  i.e. 
parallel  to  the  water  surface,  if  its  curvature  is  neglected. 

Vhen  using  the  velocity  potential  (iTa)  in  the  Bernoulli  equation  it 
must  he  rememhered  that  the  coordinates  of  cp  refer  to  a  frame  ■vdiich  moves 
with  Ax'  and  Ay*  vith  respect  to  the  rotating  frame.  Since  Bcp/dt  must  he 
referred  to  a  frame  idiich  is  at  rest  with  respect  to  the  rotating  frame, 
appropriate  convective  terms  must  he  used  to  account  for  the  motion  of 
these  frames.  The  Bemouilli  eqTiation  then  takes  the  form 

2 

_  £  «  .  ^  +  (yb)  /2  +  Ax*  +  Ay* 
p  ot  3x  3y 


(Ax*  sin  0  -  Ay*  cos  0)  -  2(b^ 

® 


2 

^(R  -  Az)^  +  2(R  -  Az)  r  cos 


0  +  r  cos' 


(l8a) 


+  ^  +  2£y  r  sin  0  +  r* 


sin® 


-  g^  r  cos  0  -  ^  r  sin  0. 

The  potential  Q  and  the  "gravitational"  caii5)onents  g^  emd  gy  do  not  refer  to 
the  gravity  of  the  earth.  Above, the  body  force  per  unit  mass  in  a  rotating 
tank  has  been  loosely  called  "gravity",  because  it  is  the  purpose  of  the 
rotating  tank  to  simulate  gravity,  g^  and  g^  are  used  in  this  sense.  Actual 
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hi^  gravity  tanks  are  built  vith  a  vertical  axis  of  rotation,  so  that  and 
gy  due  to  the  true  gravity  are  zero.  The  purpose  of  introducing  i2vas  to 
permit  an  Sitemative  discussion  of  the  effect  of  rotation.  Thus,  in  the 
above  and  the  following  relations  either  »  or  g^^  and  gy  must  be  zero. 

Insertion  of  (l8a)  into  (l6)  yields  the  following  simultaneous 
differential  equations  for  the  bubble  migration  in  a  rotating  tank  iftxen  Ax 
is  replaced  by  -Az: 


Ay”  +  3  ^  Ay*  -  2u)Az*  +  2u)^Ay  +  2gy  ■  0 


A* 


Az”  +  3  X"  *  2u)Ay*  -  2cb‘^(R  -  Az)  -  2g  «  0 


(l8b) 


2 

The  terms  with  2iu  comprise  the  Coriolis  migration,  those  with  co  the  migration 
caused  by  the  centrifugal  acceleration,  then  g^  ■  gy  ■  0. 


Setting  (I)  ■  0,  one  may  use  and  gy  in  connection  with  the  relationship 
for  the  acceleration  in  rotating  systems  (4)  to  check  this  result.  Assume 
for  a  moment  that  g^  and  gy  are  the  components  of  gravity,  i.e.  assume  a 
resting  system  with  the  x-axis  not  exactly  vertical.  Both  ^  and  slant 
downward  and  migration  by  buoyancy  is  vertically  upward.  Thus,  Ax’  and  Ay* 
are  in  the  opposite  direction  of  g^^.  and  gy  respectively.  This  is  in  agreement 
with  (l8b)  since  Ax  ■  -  Az.  For  the  migration  in  a  rotating  system  g;  and 
in  (l8b)  cannot  be  directly  replaced  by  the  centrifugal  and  Coriolis  accelera¬ 
tions  (4).  Only  the  terms  xa?  and  ym?  establish  a  "gravitational  field"  in 
the  rotating  system  idiich  causes  buoyancy.  The  Coriolis  terms  do  not  produce 
such  effects,  but  are  merely  a  device  to  describe  the  motion  of  a  body  in  a 
rotating  system  enploying  the  concept  that  the  body  bdiaves  as  if  it  were 
subjected  to  additional  forces,  the  Coriolis  forces.  To  obtain  the  mass  of 
the  body  idiich  is  of  our  interest,  namely  the  piilsating  bubble,  we  write  (l8b) 
in  the  form 


dt 


4tt 

3 


P 


4tt  A^ 
3 


0  . 


(l8c) 


This  equation  can  be  interpreted  as  an  equation  of  motion  of  a  body  ifliich  has 
the  mass  1  «  a3  and  is  subjected  to  the  force  of  buoyancy  „  4tt  a3  . 

2  3  ®y  3 

This  reflects  the  well  known  situation  that  the  virtual  mass  of  a  spherical 
bubble  in  translatory  motion  is  one-half  of  the  mass  of  the  water  displaced 
by  it.  The  Coriolis  force  ■viiich  causes  the  side  motion  of  the  bubble  is  equal 
2u)A  z*TT  1/2  p  Wilder  the  influence  of  this  term  the  bubble  moves 
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in  the  positive  y-direction.  A  direct  replacement  of  by  y*‘  according  to 
(4)  vould  have  resulted  in  a  Coriolis  migration  in  the  negative  y-direction 
and  it  vould  be  twice  as  strong.  The  above  consideration  of  the  Coriolis 
force  yields  the  same  result  as  that  obtained  by  means  of  the  e^roximate 
function 


V.  GRAVITY  MIGRATION 


As  a  first  step,  we  calculate  the  upward  migration  in  a 
resting  system  due  to  a  constant  acceleration  g.  For  this  case, 
differential  equation  (l8b)  has  the  solution 


AZ  »  2  g 


the 


(19) 


■which  has  been  fovind  independently  by  Herring  W/  and  Zoller  8/.  To  evaluate 
this  expression,  we  need  a  relationship  for  the  bubble  radius  A  as  a  function 
of  time.  The  strict  relationship  is  exceedingly  complicated,  but  according 
to  a  proposal  of  Zoller  a  satisfactory  approximation  can  be  obtained  by 

A^(t)  =  (l  +  c  -  c  (~  -  l)^)  =  A^3  c  ^ 


Here,  is  the  minimum  bubble  radius,  T,  the  period  of  pulsation,  t,  the 
reduced  time,  and 


1 


2: 

or 


1  +  c” 


Since  the  maximum  radius  is  coDmonly  substantially  larger  than  the 
minimum  radius  A^,  c  »  1  and  or  ~  At  every  place,  where  it  is  permissible, 
we  "vd.!!  neglect  c“l  in  comparison  with  xmity,  in  particular  we  can  set 
1  +  a  ~  2,  but  a  -  1  ~  l/2  c* However,  these  approximations  cannot  be  made  a 
priori,  and  each  case  must  be  considered  separately  before  the  siniplifica- 
tion  is  introduced.  In  order  to  avoid  the  presentation  of  overly  complicated 
equations  only  the  final  simplified  results  are  given  below. 

Integrating  (19)  once  we  obtain  for  the  velocity  of  migration 

AZ-  =  g  T  /  a3  dT  =  ^  (t  +  )  ,  (21) 

aP  J_^  3  Of  -  T 

•where 

a  B  A  /A  . 
m 
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The  migration  is  foxmd  to  be 
AZ(t)  = 


T  T 

I  tJ. 

-1  -1 


6  L  2 


a^  dT 


(22) 


+  2  In 


This  relation  holds  for  first  cycle  of  the  bubble  pulsation;  ve  denote 
the  bubble  period  by  T]^  and  the  value  of  a  appropriate  to  this  cycle  by  ^ 
At  the  bubble  maximum  t  »  0,  thus  the  migration  from  the  moment  of  the 
explosion  up  to  the  first  bubble  maximum  is 


AZ. 


g  Tn 


2  In-i - -  -  i  - 


(23) 


Here, 


1 


can  be  neglected  in  conrparison  with  unity  and  we  obtain  with  a  ® 
^2  2  1 
g  T^  f  g  T-i 

AZ,„  = [2  In  2-4  I  =  . .  1.5816.  (24) 


T 


The  term  gT^/8  is  the  distance  a  free  body  would  have  moved  xinder  gravity  at 
the  time  T/2^  starting  from  rest.  Thus,  the  average  acceleration  of  the 
migrating  bubble  is  I.5816  times  that  of  gravity.  A  non-p^llsating  gas  filled 
balloon  rises  in  the  water  with  an  acceleration  of  2g,  as  shown  in  equation 
(18b).  Thus,  the  expanding  bubble  rises  more  slowly  than  the  stationary 
spherical  bubble. 

Since  the  bubble  remains  spherical  during  the  expansion,  (24)  shoiild  be 
a  good  approximation  for  the  migration  during  this  phase. 


The  migration  from  the  moment  of  the  explosion  up  to  the  first  bubble 
minimum  is 


AZ, 


Im 


AZ 


g  T, 


In  4 


(25) 


This  relation  affords  a  possibility  to  study  the  dependence  on  migration  of 

Cl  or  It  is  known  that  a  bubble  contracts  to  a  smaller  size  vhexi 

under  the  influence  oT  gravity.  Thus,  A  /A  must  d^end  on  the  Froude  number 

M  m 


(26) 
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% 


To  evaluate  F,  we  use  the  following  interrelationship  between  T  and  A 


M 


(27) 


Here  t  is  the  dimensionless  period  of  the  classic  bubble  theory  and  a^,  the 
dimensionless  maximum  radius,  ll/.  is  the  total  hydrostatic  head  at  the 
^oint  of  explosion,  namely  depth  +  3^  ft  (fresh  water).  The  numerical  value 
t/a^  »  1. 6  used  in  (2?)  is  appropriate  for  a  variety  of  explosives  and  firing 
conations.  The  Froude  nvimber  then  becomes 


and  (25) 


Z 

F  .  -1 - 2 

3.84 

^^1  3.84  ^  ,  ,  ,  c 

—  ®  ^  1“  4c^/il. 5F 


Z 


(28) 


(29) 


Commonly,  bubble  migration  is  calculated  by  the  semi-empirical  formula 


1/2  , 

AZ  -  /Z^  . 


(30) 


The  constant  C  depends  on  the  property  of  the  explosive,  e.g.  for  TUT  C  a  80. 
The  equivalent  alternative  form 


.3/2 

a  0.233  F 


(31) 


is  more  general  and  holds  for  most  e3q)losives,  7/.  Comparison  of  (29)  and 
(31)  yields 

,  A,  UP 

Xn  -  1.37  Q  J3  J  .2.S8F  .  (32) 

o 


O'  =1.09, 
~  20- 


5,  the  above  relation  yields  c_  »  5.4  and  A^^/A  «  1.86.  Then, 

For  a  non-migrating  bubble  (-vdilch  is  not  covered  by  (3l))j 
Since  the  maximum  radius  is  unaffected,  this  example  illustrates 
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the  drastic  change  of  the  minimum  bubble  radius  caused  by  migration.  It  is 
also  seen  that  the  approximation  a  »  1  is  barely  acceptable  for  the  above 
conditions. 


The  rate  of  rise  of  the  bubble  in  the  second  cycle  is 

^  "^2  r  ^ 

az-2  (t)  -  -Y—  < 


2 

a  +  T 

o  J 


and  the  migration 


AZg  (t)  ^ 


-  1 


+  4  In 


«2+  1 

a  -  T 

2 


+  2  In 


a  +  T 
2 

“2  “  1 


Here,  T  is  the  period  of  the  second  cycle  and  ot  g  the  appropriate 
value  of  Qi  for  this  cycle. 


At  the  end  of  the  first  cycle  and  the  beginning  of  the  second  cycle 
the  rates  of  rise  and  AZg  (t  »  -  l)  coincide  only  if  =  Tg  and 

a  Og,  For  undam5>ed  pulsations  of  non-migrating  bubbles  these  magnitudes 
are  equal.  Energy  disipation  near  the  bubble  minimum  as  well  as  upward 
migration  change  these  magnitudes,  so  that  according  to  our  analysis  the 
rate  of  migration  changes  discontinuously  at  the  bubble  minimum.  It  is 
possible  to  remove  this  discontinuity  by  a  proper  choice  of  the  integration 
constant:  2 

AZ-  M  =  _L_  ..  (33a) 

2  3  a  -  T  -  1  T 

2.  1  2  J 


2 

T  +  — - 

a  - 

L  p 


\"2-  1 


(33a) 


AZg  (t) 


.  gii-f  -I-.  ,  — 

2  I  TgJ  “2  -  t 


+  2  — - ~^ln— ^ - 

oef-  1  2  a  -  1 

1  2 


(34a) 
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% 


For  the  piirposes  of  this  study  the  simpler  forms  given  above  are  sirfficiently 
accxirate  and  the  total  migration  up  to  the  second  bubble  maximum  becomes 


g  T 


A,Z  +  AZ 
2M  1 


4  In  2  - 


1  +  2  In  2c 

2  2 


g  T 


2  In  4c, 


and  up  to  the  second  bubble  minimuml 


AZ  +  AZ 
2  1 


In  4c  + 

2 


2  In  4c^ 


(35) 


(36) 


If  the  period  T  and  the  coefficient  c  were  the  same  in  the  first  and  second 
cycle,  the  second  cycle  migration  would  be  three  times  that  of  the  first 
cycle. 


VI.  SIDEWARD  MIGRATIOH  IN  A  ROTATING  TANK 

We  designate  the  side^Tard  coordinate  in  a  rotating  tank  by  y  and  retain 
z  for  the  upward  coordinate.  The  approximate  equation  for  the  bubble  motion 
is  according  to  the  considerations  of  Part  IV  of  this  report 


a*  2 

^y*  3'^  ^y*  -  2  00 Az*  +  2  Ayto  r:  0  .  (37) 

This  equation  can  be  readily  evaluated  following  the  previous  lines  if  the 
last  term  is  neglected  and  Az*  is  assmed  to  result  from  a  constant 
acceleration  Ikir,  yhere  R  is  the  distance  between  the  point  of  explosion 
and  the  center  of  the  rotation.  An  estimte  of  the  error  of  this  approxi¬ 
mation  ■vrill  follow. 


The  rate  of  sldeirard  motion  is, 

4  0)"  R  r  r  .3 

A^  dt  dt 


,3 


2 

=  oj-T.  T^ 


(38) 


■(diere  is  the  integration  constant  of  the  first  integration  and  K2>  that 
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of  the  second  integration.  is  readily  found  to  be 


# 


K 

1 


=  2/3 


First  cycle 


K 


=  2 


Second  cycle  , 


(38a) 


The  integration  constant  Kg  is  determined  so  that,  at  the  beginning  of  the 
second  cycle,  Ay  has  the  same  value  as  at  the  end  of  the  first  one,  if  T 
and  a  are  eq.ual.  (The  more  elaborate  treatment  shovm  before,  (33a)  and 
(34a),  does  not  appear  to  be  worth\diile  here.)  One  obtains 


K. 


K 


0 


3  ’  3  “  1 


- 1) 


2 

3 


First  cycle 


Second  cycle 


(38b) 


The  sidevard  migration  becomes 


,.:5  3 


Ayi  (t) 


(t) 


R  T  /  3 

J:  ^  T  •»-  5  +  2 

36 


12 


Ot  +  1\ 

Y  331  )  j  First  cycle 


od3  R  T^ 


(39) 


3 


2  T-  +  1  _ 

2  V  36  12 


oe  +  1 

8  ^  2  ^ 

3  -  T 


Ot  ^  r 

+  In  — ^ - r 

3  Of  _  -  1 


Second  cycle 


For  the  sideward  migration  within  the  first  and  second  cycle  respectively, 
we  find 


3  3 

0:1  R  Tt 


to  Tj^  Az^ 


4c,  - 


3  ^  ‘^''1  9 


(40) 


4 


^  3 

J  T->  m 


R  T 


1  J 


t 
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t 


t 


T.  r, 

03  R 


—  CD  T  t-Z 
3  2  2 


10  T  ,  7 

— —  In  4c^  —  -f 

3  29 


I  OD^  R 

3S  2 


(41) 


To  estimate  the  error  of  neglecting  Ay  in  (37)^  we  find  for  the  instant  of  the 
first  bubble  minimum  2 

^  203  R  T, 


aj3  Az^  +  2(d  Ay^ 


and  for  the  second  bubble  minimum 

2cd'^  R  Tg 

aD  Az*  +  ao^y  =  - - 

d  d  J 


o  2 

1  +  4c.  +  CD  T  (in  4c  -  7/*5) 
1  1  1 


(42) 


2  +  Scg  +  <I)^T2(5  In  4c^  -  7/6) 


(43) 


In  both  ca-ses  the  last  term  in  the  braces  is  negligibly  small  compared  with 
the  terms  4c]_  or  8c2,  not  only  becaoise  c  »1,  as  generally  assumed  in  our 
calculations,  but  also  because  cd%?^  is  a  small  magnitude  in  all  practical 
cases.  Since  the  last  teirnis  represent  the  magnitude  neglected  in  (37)/  the 
approximation  appea,rs  to  be  permissible.  Tlie  same  conclusion  is  reached, 
if  this  comparison  is  made  for  other  moments  of  time.  The  justification 

"  p 

of  the  second  approximation  i^iich  assumes  a  constant  acceleration  a>  R  in 
the  calculation  of  Ay  v/ill  becoEie  apparent  in  Part  VII,  iSiere  it  will  be 
s’noim  that  the  change  of  d  z  due  to  the  charge  of  acceleration  is  exceedingly 
small.  Tlie  same  result  would  be  obtained  for  Ay. 


VII.  THE  DISTORTION  OF  THE  UH-JARD  MGRATION  IN  A  ROTATIiTG 

TEST  T.'U'Hi 


According  to  the  discussion  in  Part  IV,  we  obtain  for  the  migration  in 
the  z-direction  (toward  the  mter  surface  and  the  center  of  rotation) 


t  t 

^Az(t)  =  2Rcd^  (  a”""  f  a’ 

O*^  >  DO' 


dt  dt  -2cdj 
.t  t  (T 


A  /  A-^Ay  dt  dt 
Az  dt  dt  . 


(44) 
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The  first  term  corresponds  to  the  gravity  migration  discussed  in  Section  V. 
The  second  term  is  a  result  of  the  Coriolis  acceleration^ and  the  third 
accounts  for  tlie  variation  of  the  centrifugal  acceleration  i-riLth  depth.  We 
denote  the  contribution  of  the  Coriolis  term  by  the  subscript  c  and  obtain 
with  (38): 


First  cycle 


k  4 
oi  R  T^ 


T  -  1 

240 


40 


T  +  1  +  ^ 


1 

Of,  -  T  _ 


Second  cycle 


-  - 


4  4 

cw  R  To 


T  -  1  3(t  -  1) 

^^0  40 


M|-^(of  -  D)  in-a-ll 

^5  3^2  ct  -  1 


ao+  1 

1  2 


The  term  4(a  -  l)/3  in  the  last  esqjression  stems  from  the  integration 
constant  Tdii^  assures  equality  of  the  rate  of  rise  at  the  end  of  the  first 
and  the  beginning  of  the  second  bubble  cycle.  This  term  is  shown  for 
completeness,  but  will  be  neglected  hence  forth.  The  Coriolis  contributions 
up  to  the  first  and  second  minimum  respectively  are 


^  T,  m  ^  2  2 

m  HT  fa  ,  ,  al  3<o  T, 
-T—  5  “=1  -  3 - 5 


4  4 

cn  R  T, 


co^  R  T^^ 


•2—  2ji 

15  2 


-2  = 


2  2 

29^0^  T^^ 


4  4 

Azg  -  CD  R  T^ 


Up  to  this  point,  the  integrations  needed  for  the  various  migration  terms 
have  been  simple.  But,  the  calculation  of  the  last  term  of  the  upirard 
migration  turned  out  to  be  more  complicated  and  involves  functions  'vAiich 
are  not  readily  available.  Therefore,  a  fixrther  simplification  is  made. 
For  use  in  (44)  Az  is  approximated  by 


2  „  m  2 

CD  R  T, 


'1  V  ^ 
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■where 


First  cycle 


=  2  In  2 


q  =  In  4ct 
2  ^ 

q  =  In  4c,  -  2  In  2 

3 


Second  cycle 

q^  =  2  In  2  +  2  In  4c2 

q  =  3  In  4c 
^2  2 

q  =ln4c  -21n2. 
3  2 


The  approximation  (47)  coincides  -with  (22)  and  (3^)  at  t  =  1,  0^  and  -1. 
Althou^  the  nature  of  the  functions  and  their  derivatives  differ,  the 
approximation  (47)  id-ll  suffice  for  the  rotigh  estimates  'vdiich  are  the 
purpose  of  this  study. 

The  migration  due  to  the  change  of  the  centrifugal  acceleration  then 
becomes  for  the  i-th  cycle 


«.2t2 


••  Az. 


iz 


k  k 
(XT  RT 


12 


[(Bo 


Az  dT  dT 


(1^8) 


^  ,  a.+  1  a  +  T 

^  +  K)  In  -3 -  +  K  In  (x  +  l) 

oj  -  T  oe  1  ^ 

i  i 


+  B  (t2.  1)  +  B  (t3  +  1)  +  B.  (x^  -  1) 


id.th 


Bq  =  2q^  -  4)/l5 


(16 

In 

2  - 

4  + 

2  In  4c^)/l5 

(16 

In 

2  - 

4  + 

22  In 

“  ^ 

c 

J" 

1 

In 

4c 

"  4 

1 

-  i 

In 

4c 

4  2 


(49a) 

First  cycle 
Second  cycle 

(J^9b) 

First  cycle 
Second  cycle 
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The  total  migration  in  the  rotating  tank  is 

DAz.  ■  Az.  +  Az.  +  Az.  .  (51) 

1  i  xc  xz 

VIII.  DISCUSSION  OP  THEORBriCAL  AND  EXPERIMEHTAI.  RESULTS 

A  fev  explosion  tests  made  with  the  NOL  centrifuge  made  the  observation 
of  the  sideward  migration  possible  £/.  Unfortunately,  ah  accurate  evalua¬ 
tion  of  this  migration  proved  to  be  difficult,  so  that  the  values  given 
below  are  rather  uncertain.  Nevertheless,  these  tests  may  serve  as  an 
illustration  of  the  veurious  migration  components. 

The  test  conditions  of  one  of  these  tests  were  as  follows: 

R  »  10  ft 

,  -1 

CD  ■  14.2  sec  . 

Fig.  4  is  a  contour-time  plot  of  the  bubble  showing  the  upward  and  sideward 
migrations.  The  ciurves  give  the  two  extremes  of  the  bubble  dimension  in 
horizontal  and  vertical  extension  as  a  functicna  of  time.  The  bubble  center, 
asstmed  to  be  halfway  between  these  extremes,  is  shown  as  a  dotted  curve. 

The  graph  shows  a  difficulty  ■viiich  is  typical  of  migrating  bubble 
problems,  namely  the  xmcertainty  of  the  time  of  the  bubble  minimum.  The 
maximum  contraction  in  vertical  direction  occurs  at  a  different  time  (4.85 
msec)  than  that  in  the  horizontal  direction  (5.25  msec).  The  first  instant 
corresponds  to  the  mcsnent  vhere  the  lower  and  upper  bubble  interface  impinge 
on  each  other.  In  pressure  records,  this  Instant  is  manifested  by  a  short, 
but  strong  pressure  pulse  originated  by  the  water  hammer.  This  "spike”  is 
superposed  on  the  bubble  pulse  originated  by  the  conpressed  gases  in  the 
bubble  and  it  occurs  before  the  pressure  maximum  of  the  bubble  pulse,  Pig.  5* 
The  latter  moment  corresponds  to  that  of  the  bubble  minimum  and  it  is  not 
necessarily  identical  with  the  moment  of  the  maximum  lateral  contraction. 

This  is  nicely  demonstrated  by  the  dotted  ctuve  r^resenting  the  upward 
migration  in  Pig.  4:  Significant  upward  migration  (i.e.  regions  ^ere  the 
rate  of  migration  is  hi^)  begins  at  t  ■  4.85  msec  (vertical  minimum  of 
bubble)  and  seems  to  end  at  t  ■  5.25  msec  (horizontal  minimum).  Since 
significant  migration  takes  place  before  as  well  as  after  the  moment  of  the 
bubble  minimum,  the  minimum  must  occur  between  t  ■  4.85  msec  and  t  «  5.2 
msec.  We  have  chosen  t  ■  5.1  msec  and  obtain: 


4 
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Observed  Bubble  Periods  Total  "Upvard"  Migration 

Tj^  ■  5*1  nisec  »  O.65  inch 

T  "  5*2  msec  XAz  +  *»  I.60  inch 

2  2 

Since  the  difference  betveen2)Azj_  and  Azj^  is  small,  ve  set  approximately 

Az^  “  0.65  inch 

Azg  “  0.95  inch. 

From  (4o)  and  (4l),  ve  find  for  the  sidevard  motion 

Ay^  a  0.03  inch 

Ay^  “  0.098  inch. 

The  smallness  is  these  magnitudes  shovs  that  a  quantitative  comparison  vith 
measured  data  is  a  doubtful  undertaking,  since  3/IOO  or  even  l/lO  of  an 
inch  are  not  readily  evaluated  from  a  hi^-speed  movie  film.  Nevertheless, 
such  an  evaluation  is  atteii5)ted  in  Pig.  6  ’Hhich  shows  an  acceptable  agreement 
between  measured  and  calculated  side  migration.  At  the  end  of  the  second 
cycle  the  calculated  side  migration  is  about  30^  higher  than  the  experimental 
one,  0.1  vrs,  0.13  inch.  This  would  be  in  line  with  other  theoretical  results 
based  on  spherical  bubble  shapes.  Because  of  possible  baseline  shifts,  the 
side  migration  at  the  end  of  the  first  cycle  is  uncertain.  However,  this 
comparison  indicates  that  our  calculations  give  a  satisfactory  estimate 
of  the  sideward  migration. 

The  Coriolis  correction  for  the  "upward"  migration  as  obtained  from 
(46)  is: 

Az^  =  -  1.5  '  10”^  Az., 

Ic  1 

Azgc  »  -  T.10“^  Az^ 

and  the  correction  because  the  variation  of  the  acceleration  is 
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0. 6 


10"^  Az, 


-  1.3  •  lO"^  Az  . 

2 


It  is  interesting  that  the  two  "errors"  have  the  same  sign.  However,  the 
error  due  to  change  of  acceleration  is  much  smaller  than  that  caused  by  the 
Coriolis  effect.  The  total  error  in  "upmrd"  migration  is  0.8  percent, 
hence  negligibly  small  for  practical  purposes  idiere  experimental  errors  of 
much  larger  magnitude  are  to  be  expected. 


Since  the  conditions  of  this  example  are  typical  of  a  high  gravity  tank 
setup,  one  is  justified  in  concluding  that  the  errors  in  the  upt'fard  migration 
of  explosion  bubbles  which  are  introduced  by  the  rotation  will  not  constitute 
an  obstacle  against  adopting  this  type  of  high  gravity  tank. 

Equations  (46)  and  (50)  show  that  the  errors  in  upward  migration  involve 
terms  ^ike  to  T  ,  cd  T  H/Az,  etc.  With  g^  *  the  first  of  these  terms 

is  g^T  /R.  Hence,  in  this  case  the  error  term  is  inversely  proportional  to 
R,  for  any  given  value  of  g^T  .  As  the  ot^er  terms  show,  the  total  error 
is  not  necessarily  a  linear  function  of  R  ,but  involves  higher  negative 
powers  of  R.  This  confirms  the  argument  that  a  slowly  rotating  centrifuge 
with  a  large  am  produces  less  distortion  of  the  bubble  migration  than  a 
small,  fast  rotating  centrifuge. 


author  thanks  ttr.  R.  S.  Price  for  the  preparation  of  a  graph  on  which 
Figure  4  on  this  report  Is  based.  The  efforts  of  Dr.  E.  Swift,  Jr.  who  has 
read  the  draft  of  this  report  are  appreciated. 
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FIG.  I  COMPONENTS  OF  THE  VECTOR  [oi't] 

(the  vector  [wCwT]]  is  not  shown) 


AXIS 
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WATER  SURFACE 


FIG.  2  TRAJECTORIES  OF  SPRAY  PARTICLES  WHICH  MOVE 
INITIALLY  TOWARD  THE  CENTER  OF  ROTATION,  AS 
SEEN  BY  A  ROTATING  OBSERVER 
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FIG.  4  BUBBLE  CONTOURS  AND  MIGRATION  IN  UPWARD  AND  SIDEWARD  DIRECTION 


NOLTR  61  -  145 


* 


UJ 

X 


h- 

.  o 

UJ  CO 

CO  o 


3  UJ 
CL  -J 
CD 
UJ  00 


LU 

X  U. 

H  O 


Ul 

o 


2 

O 


UJ  H- 
CL  < 
<  QC 
X  3 
CO  Q 


10 

O 

u_ 


28 


NOLTR  61  -  145 


29 


FIG.  6  SIDEWARD  MIGRATION 
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